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ABSTRACT 

We have recently obtained the exact analytic solutions of the relativistic equations relating the radial and 
time coordinate of a relativistic thin uniform shell expanding in the interstellar medium in the fully radiative 
and fully adiabatic regimes. We here re-examine the validity of the constant-index power-law relations between 
the Lorentz gamma factor and its radial coordinate, usually adopted in the current Gamma-Ray Burst (GRB) 
literature on the grounds of an "ultrarelativistic" approximation. Such expressions are found to be mathemati- 
cally correct but only approximately valid in a very limited range of the physical and astrophysical parameters 
and in an asymptotic regime which is reached only for a very short time, if any, and are shown to be not 
applicable to GRBs. 

Subject headings: gamma rays: bursts — ISM: kinematics and dynamics — relativity 



1. INTRODUCTION 

The discovery of afterglows (ICosta et alJ Il997j) has of- 
fered a very powerful tool for reaching an understanding of 
Gamma-Ray Bursts (GRBs). The consensus has been reached 
that such an afterglow originates from a relativistic thin shell 
of baryonic matter propagating in the InterStellar Medium 
(ISM) and that its description can be obtained from the rel- 
ativistic conserva t ion la ws of energy and momentum. In 
iBianco & Ruffinil J2005I) we reported the exact analytic so- 
lutions of the corresponding equations, respectively under 
fully radiative and fully adiabatic conditions, giving in both 
cases explicit relations between the laboratory time and the 
radial coordinate of the shell. We here compare and contrast 
our results with the simple constant-index power-law rela- 
tion between the Lorentz gamma factor and the radial coor- 
dinate of the shell generally adopted in the current literature, 
obtained usin g a so -called "ultrarelativist ic" approximation 
(see e.g. I^ll997i fl99l IWaxmanlll997t iRees & Meszaros! 
" 1998HGranot et al ll999| lPanaitescu & Meszaros 1998 Jl999t 
Chian g & Permed 1 19991: iPiranl 1 19991 iGruzinov & Waxmanl 
1999; van Paradiis et al. 2000; Meszaros 2002, and references 



therein). We show that such an approximation only holds in 
a very limited range of the physical and astrophysical param- 
eters and in an asymptotic regime which is reached only for 
a very short time, if any. We demonstrate that this constant- 
index power-law cannot be used for modeling GRBs. Illustra- 
tive examples are given for the source GRB 991216. 

2. THE AFTERGLOW ANALYTIC SOLUTIONS 

The fulfillment of the energy and momentum conservation 
for the equations of motion of the relativistic baryonic matter 
shell in the labora tory reference frame leads to the following 
equations (see e.g. Piran 1999; Ruf fini et all2003bl and refer- 
ences therein): 

dE int = (j-l)dM ism c 2 , (la) 
d 1 = -[( 1 2 -\)/M] dM ismi (lb) 
dM=[(l-e)/c 2 ]dE int + dM lsm , (lc) 

1 ICRA - International Center for Relativistic Astrophysics. 

2 Dipartimento di Fisica, Universita di Roma "La Sapienza", Piazzale Aldo 
Moro 5, 1-00185 Roma, Italy. E-mail: bianco@icra.it, ruffini@icra.it 
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n r 2 dr, 



(Id) 



where M is the shell mass-energy, n; sm is the ISM number 
density, m p is the proton mass, e is the emitted fraction of the 
energy developed in the collision with the ISM and M; sm is 
the amount of ISM mass swept up by the shell within the ra- 
dius r. Mi sm = (47r/3)m ; ,«i sm (r 3 -r 3 ), where r is the starting 
radius of the baryonic matter shell. In general, an additional 
equation is needed in order to express the dependence of e 
on the radial coordinate. In the following, e is assumed to be 
constant and such an approximation appears to be correct in 
the GRB context. 

Consensus has also been reached on a simple integration 
of the equations of mo tion Eqs.© in the f ully radiative case 
(s= 1 ) fsee lPiradl999l:lRuffini etalJ2003btlBianco & Ruffinil 
2005)), leading to: 



7 = 



1 + (Mjsm/Mfi) (1 +7o') [1 + \ (Mism/Mj,)] 

Y" 1 + (M ism /M B ) (1 +7o') [1 + 2 (M ism /M B )] 



(2a) 



where M B and j a are the initial values respectively of the 
mass and of the Lorentz gamma factor of the baryonic 
sh ell. The above Eqs.([T}.(|2al differ from the ones derived 
by Bl andford & McKeef 7l976) in a different framework, of- 
ten quoted in the current literature within the present con- 
text. Correspondingly, in the fully adiabatic c ase (e = 0), 
Eqs.lJD have the following analytic solution (see Piran 1999; 
Bianco & Ruffinil2005h 1: 



2 _ 7g + 2 7o (Mism/Mfl) + (M ism /M B ) 
1 + 2 7o (M ism /M B ) + (M ism /M B ) 1 



(2b) 



In IBianco & Ruffinil J200l we have explicitly integrated 
the differential equations for r(t) in Eqs.iyJ, recalling that 
7 -2 = 1 - [dr/(cdt)] 2 , where t is the time in the laboratory 
reference frame. We have then obtained new explicit ana- 
lytic expressions for the equations of motion of the relativis- 
tic shell which are essential for explicitly obtaining the ana- 
lytic expressions of the equitemporal surfaces (EQTSs) in the 
full y radiative and in the fully ad iabatic cases, respectively 
fsee lBianco & Ruffina200ll2005h . 

3. APPROXIMATIONS ADOPTED IN THE CURRENT LITERATURE 
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We turn now to the comparison of the exact solutions given 
in Eqs.0 with the approxima tions used in t he current litera- 
ture. Following Blandford & McKee ( 1976), a so-called "ul- 
trarelativistic" approximation 70 3> 7 3> 1 has been widely 
adopted by many authors to solve Eqs.{l} (see references in 
section^. This leads to simple constant-index power-law re- 
lations: 

7 r~ a , (3) 

with a = 3 in the fully radiative case and a = 3/2 in the fully 
adiabatic case. This simple relation is in contrast with the 
complexity of Eqs.(|2|i. 

We address now the issue of establishing the domain of ap- 
plicability of the simplified Eq.(|3} used in the current litera- 
ture both in the fully radiative and adiabatic cases. 

3.1. The fully radiative case 
We first consider the fully radiative case. If we assume: 



l/(7„ + l)«M lsm /Mj«7o/(7o+l)< 1- 



(4) 



we have that in the numerator of Ea.d2al the linear term in 
Mis m /MB is negligible with respect to 1 and the quadratic term 
is a fortiori negligible, while in the denominator the linear 
term in M^/Mb is the leading one. Ea.d2al> then becomes: 



[7o/(7o + l)]M B /M ism . 



(5) 



If we multiply the terms of Eq.® by (70 + l)/7o, we obtain 
l/7o < (M ism /M B )[(7o + l)/7o] < 1, which is equivalent to 
7o > [7o/(7o+ l)](M B /M ism ) > 1, or, using Eq.(0, to: 



7o > 7 > 1 . 



(6) 



which is indeed the inequality adopted in the "ultrarelativis- 
tic" approximation in the current literature. If we further as- 
sume r 3 ^> r\, Eq.Q can be further approximated by a simple 
constant-index power-law as in Eq.(|3j: 



7 ~ [ 7o /( 7o +l)]M B /[(4/3) 



7rn ism m„r 3 l oc r 3 . (7) 



We turn now to the range of applicability of these approx- 
imations, consistently with the inequalities given in Eq.@. 
It then becomes manifest that these inequalities can only be 
enforced in a finite range of M; s111 /Mb. The lower limit (LL) 
and the upper limit (UL) of such range can be conservatively 
estimated: 

(tr) LL = lo2 ^r> (t) UL = 10 ~ 2 ^- (8a) 

The allowed range of variability, if it exists, is then given by: 

{t) UL -{t) LL = 10 ~ 22 ^ >0 - (8b) 

A necessary condition for the applicability of the above ap- 
proximations is therefore: 



7o > 10 4 



(9) 



It is important to emphasize that Eq.© is only a necessary 
condition for the applicability of the approximate Eq.Q, but 
it is not sufficient: Eq.0 in fact can be applied only in the 
very limited range of r values whose upper and lower lim- 
its are given in Ea.d8ali. See for explicit examples section 
below. 



3.2. The adiabatic case 

We now turn to the adiabatic case. If we assume: 

l/(2 7o )<M ism /M B <7o/2, (10) 

we have that in the numerator of Ea. J2b> all terms are negli- 
gible with respect to 7^, while in the denominator the leading 
term is the linear one in Mj sm /M#. Eq. (l2b> then becomes: 



7~^( 7o /2)M B /Mi 5 



(11) 



If we multiply the terms of Ea. (llO> by 2/70, we obtain 
I/7 2 < (2/7o)(M ism /M B ) < 1, which is equivalent to 7 2 > 
(7o/2)(M B /M ism ) > 1, or, using Eq.O, to: 



7 2 > 7 2 > 1 . 



(12) 



If we now further assume r 5 r 3 , Eq.l^J can be further 
approximated by a simple constant-index power-law as in 
Eq.(3): 



7 ~^/( 7 o/2)M fl / [(4/3) 



irn ism m p r 3 ] oc r 3/2 . 



(13) 



We turn now to the range of applicability of these approx- 
imations, consistently with the inequalities given in Eci.lllOL 
It then becomes manifest that these inequalities can only be 
enforced in a finite range of M\ sm /M B - The lower limit (LL) 
and the upper limit (UL) of such range can be conservatively 
estimated: 

(*L- 10r **- «*> 

The allowed range of variability, if it exists, is then given by: 

(14b) 



=10- 2 ^ >0. 

Mb J,, 270 



A necessary condition for the applicability of the above ap- 
proximations is therefore: 

7o>10 2 . (15) 

Again, it is important to emphasize that Eq.J15i is only a 
necessary condition for the applicability of the approximate 
Eq.jl3L but it is not sufficient: Eq.jl3> in fact can be applied 
only in the very limited ran ge of r values whose upper and 
lower limits are given in Eq.( ll4at . See for explicit examples 
section|4]below. 

4. A SPECIFIC EXAMPLE 

Having obtained the analytic expression of the Lorentz 
gamma factor for the fully radiative case in Eq.d2ali, we illus- 
trate in Fig. ^the corresponding gamma factor as a function 
of t he radial coordinate in the afterglow era for GRB 991216 
(see Ruffini et al. 2003b, and references therein). We have 
also represented the corresponding solution which can be ob- 
tained in the adiabatic case, using Eq.( l2bt . starting from the 
same initial conditions. It is clear that in both cases there is 
not a simple power-law relation like Eq.(|3} with a constant 
index a. We can at most define an "instantaneous" value a e /f 
for an "effective" power-law behavior: 



2 eff : 



c/ln7 
(ilnr 



(16) 



Such an "effective" power-law index of the exact solution 
smoothly varies from to a maximum value which is always 
smaller than 3 or 3/2, in the fully radiative and adiabatic cases 
respectively, and finally decreases back to (see Fig. Q. We 
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FIG. 1 . — In the upper panel, the analytic behavior of the Lorentz 7 factor 
during the afterglow era is plotted versus the radial coordinate of the expand- 
ing thin baryonic shell in the fully radiative case of GRB 991216 (solid red 
line) and in the adiabatic case starting from the same initial conditions (dot- 
ted blue line). In the lower panel are plot ted t he corresponding values of the 
"effective" power-law index a e ff (see Eg. 1161 ). which is clearly not constant, 
is highly varying and systematically lower than the constant values 3 and 3/2 
purported in the current literature (horizontal dotted black lines). 



see in particular, from Fig. [2 how in the fully radiative case 
the power-law index is consistently smaller than 3, and in the 
adiabatic case a e ff = 3/2 is approached only for a small inter- 
val of the radial coordinate corresponding to the latest parts 
of the afterglow with a Lorentz gamma factor of the order of 
10. In this case of GRB 991216 we have, in fact, 7o = 310.13 
and neither Eq.(|6j nor Eq. (I12> can be satisfied for any value 
of r. Therefore, neither in the fully radiative nor in the adia- 
batic case the constant-index power-law expression in Eq.Q 
can be applied. 

For clarity, we have integrated in Fig. [2] an ideal GRB 
afterglow with the initial conditions as GRB 991216 for se- 
lected higher values of the initial Lorentz gamma factor: 
7o = 10 3 , 10 5 , 10 7 , 10 9 . For 7o = 10 3 , we then see that, again, 
in the fully radiative condition a e /f = 3 is never reached and 
in the adiabatic case a e /f ~ 3/2 only in the region where 
10 < 7 < 50. Similarly, for 70 = 10 5 , in the fully radiative 
case a e ff ~ 3 is only reached around the point 7 = 10 2 , and 
in the adiabatic case a e ff ~ 3/2 for 10 < 7 < 10 2 , although 
the non-power-law behavior still remains in the early and lat- 
est afterglow phases corresponding to the 7 = j a and 7 — ► 1 
regimes. The same conclusion can be reached for the remain- 
ing cases 70 = 10 7 and 7 Q = 10 9 . 

We like to emphasize that the early part of the afterglow, 
where 7 = 70, which cannot be described by the constant- 
index power-law approximation, do indeed corresponds to the 
rising part of the afterglow bolometric luminosity and to its 
peak, which is reached as soon as the Lo rentz gamma factor 
starts to decrease. We have shown (see e.g. Ruffini et al. 2001 , 
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2003b] 120051 and references therein) how the correct identi- 
fications of the raising part and the peak of the afterglow are 
indeed crucial for the explanation of the observed "prompt 
radiation". Similarly, the power-law cannot be applied dur- 
ing the entire approach to the newtonian regime, which cor- 
responds to some of the actual observations occurring in the 
latest afterglow phases. 

5. CONCLUSIONS 

It is well known that scaling laws and constant-index 
power-law expressions are obtainable only in the asymptotic 
case of ultrarelativistic regimes and in the Newtonian limit, 
while in the f ully relativisti c regime the scaling laws break 
down (see e.g. Rul'linil 1 97 31) . This circumstance is more sub- 
tle in GRB afterglows: 1) the ultrarelativistic approximation 
is only a necessary condition, but not a sufficient one, for the 
existence of scaling laws; 2) such a necessary condition im- 
plies values of the initial Lorentz gamma factor 70 outside the 
range currently observed in G RB s ources. 

We have shown in section 13.11 that, in the fully radiative 
case, the necessary ultrarelativistic condition for obtaining 
the appearance of scaling laws is 70 > 10 4 . We recall that 
the 70 values deduced typically for GRBs are of the order of 
7o ~ 10 2 (see e.g. GRB 030329, GRB 02 0322, GRB 991216 
GRB 98 0519, GRB 980425, GRB 970228 Ruffini et al. 2002, 
2003b a, and references therein). Thus this necessary condi- 
tion is never fulfilled in GRBs. 

It would appear from section l3~2l that the necessary ultra- 
relativistic condition for obtaining the appearance of scaling 
laws is less severe in the adiabatic case: 7 > 10 2 . How- 
ever, this condition is not sufficient for the applicability of the 
constant-index power-law approximation to the entire after- 
glow, as clearly shown in Fig. ^ The regime a e ff = 3/2 is 
in fact approached only asymptotically and anyway in a very 
limited region. 

In the current literature (see references in section[0, a sys- 
tematic use of the constant-index power-law approximation 
for the Lorentz gamma factor has been made. The 7 = 70 
regime has been generally neglected or erroneously matched 
to the constant-index power-law approximation hampering 
the un derstanding of the observe d "prompt radiation" (see e.g. 
Ruffin i et al J20 0lL l20^ l2005l). 

We expect that the data from Swift will soon add observa- 
tional evidence to the validity of this theoretical treatment. 
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FIG. 2. — In these four diagrams we reproduce the same quantities plotted in Fig. Qfor four higher values of 70 . The upper (lower) left diagram corresponds 
to 70 = 10 3 (70 = 10 5 ). The upper (lower) right diagram corresponds to 70 = 10 7 (70 = 10 9 ). It is manifest how asymptotically, by increasing the value of 70, 
the values a = 3 and a = 3/2 (horizontal black dotted lines) are reached, but only in a limited range of the radial co-ordinate and anyway for values of 70 much 
larger than the ones actually observed in GRBs. 
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